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Abstract 

The fluctuations of the flux radiated by an evaporating black hole will be discussed. 
Two approaches to this problem will be adopted. In the first, the squared flux operator 
is defined by normal ordering. In this case, both the mean flux and the mean squared 
flux are well defined local quantites. It is shown that the flux undergoes large fluc- 
tuations on a time scale of the order of the black hole's mass. Thus the semiclassical 
theory of gravity, in which a classical gravitational field is coupled to the expectation 
value of the stress tensor, breaks down below this time scale. In the second approach, 
one does not attempt to give meaning to the squared flux as a local quantity, but only 
as a time-averaged quantity. In both approaches, the mean squared mass minus the 
square of the mean mass grows linearly in time, but four times as fast in the second 
approach as in the first. 

1 Introduction 

One of the most remarkable theoretical discoveries of recent decades was that of black hole 
evaporation by Hawking [jl], |2| in 1974. This discovery demonstrated that a black hole emits 
a (filtered) thermal spectrum of radiation, and will eventually cease to exist as a classical 
black hole. However, Hawking's derivation of this effect and most of the subsequent papers 
on the topic have dealt only with the mean flux of radiation emitted by the black hole. 
There should in addition be fluctuations of the flux, which will be the topic of this paper. 

In order to discuss the fluctuations in the components of the stress tensor, it is necessary 
to be able to define the expectation value of the square of a stress tensor component. As will 
be discussed in the following section, this is a more difficult problem than the definition of 
the expectation value of the stress tensor operator itself. In flat spacetime, one possibility is 
to normal order the product of stress tensor operators. This gives a meaning to quantities 
such as the square of the energy density or pressure at a spacetime point and hence to 
a local measure of the fluctuations in these quantites. It will be shown in Sect. || that 
normal ordering here involves dropping both a divergent vacuum term and a divergent, 
state-dependent, cross term. The normal ordering approach was adopted in Refs. 0, 0, ||, 
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and was generalized to static curved spacetimes by Phillips and Hu 0. It leads to the 
correct classical limit in that the expectation value of a squared stress tensor component is 
the square of its expectation value in the limit that the quantum state is a coherent state. 
Furthermore, it leads to the prediction that quantum states which exhibit negative energy 
densities have large energy density fluctuations ||. 

Black hole evaporation necessarily involves negative energy density in that a flux of 
negative energy going down the horizon is needed to account for the mass loss by the black 
hole. Thus the results of Ref. []|] lead us to suspect that there must be large instantaneous 
fluctuations in the Hawking flux. Such large fluctuations can also be inferred on statistical 
physics grounds, as will be discussed in Sect. In the present paper, we will restrict our 
attention to the Hawking flux in the asymptotic region far from the event horizon where 
spacetime is approximately flat. In Sect. ||, this will be done using the normal ordering 
prescription. In Sect. |3.2| , a formula will be derived for the squared flux radiated by a 
moving mirror in two dimensional spacetime. This can be used to discuss the fluctuations 
in the flux enitted by a black hole in two dimensions. In Sect. 3.4, the corresponding 
analysis will be carried out for a Schwarzschild black hole in four dimensional spacetime. 
In both cases it will be shown that there are fractional flux fluctuations of order unity on 
time scales of the order of the mass of the black hole. 

However, there are alternatives to the normal ordering method which involve space or 
time averages of the stress tensor. One such alternative was used by Barton to study the 
fluctuations of the Casimir force. In Sect. || the possibility of averaging the flux over time 
in a two dimensional model will be examined. Using an integration by parts procedure, it 
is possible to define time integrals of the state-dependent cross term. It will be shown that 
the mean squared mass of the black hole undergoes a random walk in that its deviation 
from the square of the mean mass grows linearly in time. This is true regardless of whether 
or not the cross term is retained, but the rate of growth due to the cross term is three times 
that obtained in the normal ordering prescription. 

The results of the paper will be summarized and discussed in Sect. |5[ 



2 Normal-ordering expansions and the cross term 

In Minkowski spacetime we renormalize the expectation value of the energy-momentum 
tensor by subtracting out the Minkowski divergence, 

{^^i/} renormalized (• T au .) 

= (T UU )-(T UU ) M . (1) 

Here () m denotes the expectation value in the Minkowski vacuum state. For quadratic 
operators, this subtraction is just normal ordering. However, this is not true for the squared 
energy-momentum tensor. In this paper we will study the massless, minimally coupled 
scalar field, for which the energy-momentum tensor is 

T nu = <P(,n.4>,v) - \guv(t>, P 4>' p = \{4>,n4>,u + <P,u<P,u) - \guu<t>, P 4> ,p - (2) 
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The flux operator T rt = —T rt has finite expectation values. The product of a pair of normal 
ordered quadratic operators, (: 0102 0304 :) can be expressed using Wick's theorem as 

: X 2 ■■ 0304 : = : 01020304 : 

+ : 0103 : (0204)m+ : 0104 : (0203)m 
+ : 0203 : (0i04)m+ : 0204 : (4>\<fa)M 

+ (0103)a/(0204)m + (0104)Af (0203)M • (3) 

The first term is fully normal-ordered, the next four are cross terms and the final two are 
pure vacuum terms. The flux two-point function for an arbitrary state can written as 

(T rt (x)T rt (x')) 
= (:T rt (x)::T rt (x'):} 

= {: T rt (x)T rt (x') :) + (T rt (x)T rt (x')) cross + {T rt (x)T rt (x')) M , (4) 
where the vacuum term 

{T rt {x)T r t(x')) M = {4>( x ),r(p(x'),r')M{4>( x ),t(p(x'),t')M (5) 

is the expectation value in the Minkowski vacuum state and 

{T rt (x)T rt {x')) cross = {: 0(x) )r 0(x').r' ■) {<t>{x) .t4>( x '),t> ) M + (■ 0(» ,t<t>(x') ,t> ■) {4>{x) .r4>{x') ,r') M 

(6) 

is a state-dependent cross term. In the coincidence limit, x' — ► x, (T Tt {x)T rt {x')) m and 
(T rt {x)T rt (x')) cr0 ss are divergent. We can see that although the off-diagonal components of 
energy-momentum tensor are finite, their square is divergent and remains so even if (T^m 
is subtracted. We have cross terms which contain state-dependent divergences. If we wish to 
give meaning to the notion of the squared flux as a local quantity, it is necessary to remove 
this divergence. One possibility is simply to remove both the vacuum and cross term and 
define (T^ t ) = (: T^ t :). In the following section, we will consider the normal ordered term 
alone, and in Sect. ^, we will return to the issue of whether there is a nontrivial physical 
content to the cross term. 

Note that the quantities we will investigate, (: T rt {x)T rt {x') :) and (: T rt (x)T r t(x') :) 
+{T r t(x)T r t(x')) cross, are distinct from the stress tensor correlation function given by 

<W = (T^(x)T pa (x')) - (T^(x))(T po -(x')) . (7) 

This latter quantity is independent of the choice of renormalization in the sense that it 
is unchanged by subtracting a c- number from T pu , but it is singular in the coincidence 
limit. This correlation function was used by Muller and Schmid || to discuss cosmological 
perturbations due to quantum fields and by Carlitz and Willey Q in the context of black 
hole evaporation. However, for the questions which we wish to pose, it seems more natural 
to examine the various terms in Eq. (||). 
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Figure 1: The worldline of a mirror moving to the left. The radiation emitted to the right 
is illustrated. 



3 Flux fluctuations using the normal-ordering scheme 

3.1 Formula for (: T r 2 :) 

Consider the massless scalar field, whose stress tensor is given by Eq. (|2|). The normal- 
ordered expectation value of a product of fluxes is shown in Appendix A to be 



(: T rt (x)T rt (x') :) = (: T rt (x) :)(: T rt {x') :) + (: 0, r (x)0, r (x') :)(: <f>, t (x)<f>,t(a/) :) 

+ (: ^r(x)Mx') ■■)(■■ Mx)M*') ■■) (8) 

and the squared flux is 

(: T 2 (x) :) = Um x ,^ x (: T rt (x)T rt {x') :) (9) 
3.2 Two dimensional moving mirror 

In flat space-time, boundaries induce vacuum energy and stress. If the boundaries move, 
then particles can be created. A single reflecting boundary (mirror) can create particles if 
it undergoes non-uniform acceleration. 

We follow the treatment of Fulling and Davies 1C, 11], and consider a massless scalar 



field in two dimensional flat spacetime with an arbitrary mirror trajectory, as illustrated in 
Fig. 1. 

The trajectory for the mirror is 

x = z(t), (10) 
where |^'(t)| < 1, and z = for t < 0. The boundary condition for the scalar field <fi is 
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<f>(t,z(t)) = Q, (11) 
and the positive frequency mode function for t > is given by 

<f, u (x) = i(4vrw)- 1 / 2 [e-^ - e -M^u-u)^ (12) 
Here u = t — x and v = t + x are null coordinates, and the parameter t u is defined by 

Tu ~ z(t u ) = u. (13) 
The phase of the reflected wave is a function of u only and is defined by 

p(u) = 2t u - u. (14) 

The phase change of the out-going mode is due to the Doppler shift at the moving mirror. It 
is not surprising to see that the moving mirror can create particles, since the mirror bound- 
ary condition changes a positive frequency mode into a mixture of positive and negative 
frequencies. 

The quantum field operator can be written as 

roo 

<f>(x,t) = du[a u (pu + alfa], (15) 
Jo 

where a u and are annihilation and creation operators, respectively. The in-vacuum state 
|0) is defined by 

ajO) =0, V uj. (16) 

We use the point-splitting method to extract out the Minkowski vacuum divergence. Let 
</>* be replaced by + e, x), and 

POO 

(<P,n^,v) = hm / <A^(r,i)0* (r,i + e)«L;. (17) 
e ^o Jo 

The derivatives of the mode functions become 



, . \ 1/2 

^-j [e~ iujv -p'{u)e~ lu> ^ u \ (18) 
ui \ 1 / 2 

— [e~ iuJV +p'(u)e- iujp ^}, (19) 
4ir ) 

uj \ 1 / 2 

- [ e ^(«+ £ ) - p'( n + e)e itJp( - u+e \ (20) 



/ \ 1/2 

0* r= ( j [ e *"(*+e) +p '( u + e ) e ^p(«+*)]. (21) 
We insert these expressions into Eq. (|l7|) , and evaluate the w-integration using 



and 
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and 



f°° e ibx dx= 1 -, 
Jo o 



uje iauJ dLu = — / e iaw dw 
o ioa Jo 



Expanding the results in a power series in e yields 



1 


"l 


/p» 


4ir 




U' 


1 




(V' 


47T 


I 




1 

47T 


' 2 


+ ( 



6 p' 

lp'" 
6 p' 



+ 0(e), 
+ 0(e), 



2p' 



1 /p" 



and 



1 

47T 



(■u — p(u)) 2 4 
2p' 



+ 



\jT_ 
6 p' 



(f — p(it)) 2 



4VP'/ 6p' 



(22) 



(23) 



(24) 
(25) 
(26) 

(27) 



Here (<fi tT <fi t t) and (0,t0,r) are finite, and (0 jr jr )and {(/) can be renormalized by 
discarding the e -2 term. This term is independent of p(it) (i.e. independent of trajectory) 
and can be recognized as the Minkowski vacuum divergence. The normal-ordered operator 
products are 



(: <f>,r<f>,t : ) 
(: <P,t4>,r ■) 



(<j>,r<j>,t), 



+ 



2TT€ 2 ' 



and 



2TT€ 2 ' 



(28) 
(29) 

(30) 



(31) 



We next substitute the above relations into Eq. (||). The squared flux for an arbitrary 
trajectory becomes 



and the flux is 



-V 

4vr/ 



<P') 



i\2 



+ 



3 (p"\ 4 \{p") 2 p'" (p 



J/f\2 



+ 



(u-p(u)) 4 16 \p' J 4 (p') 3 12(p') 2 



(32) 



i _ l£_ 
4Vp'/ 6 p' 



(33) 
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3.2.1 A trajectory which produce a thermal spectrum of particles 



Fulling and Davies [10, |ll]] have discussed a particular mirror trajectory which produces a 
steady thermal flux of particles at late times, and which models a two-dimensional evapo- 
rating black hole. Carlitz and Willey M have shown that the correlation function C defined 
in Eq. (^) in this case is just that for a thermal state. The trajectory has the asymptotic 
form 



z(t) ~ -t-Ae~ 2Kt+B , t 
where A, B, and k are positive constants. Here 



oo, 



p(u) = 2t u -u~B- A- k{ - u+b ^ , u ^ oo. 
Substituting this form Eq. (^|) gives the normal-ordered squared flux 



(34) 



(35) 



(■■ T&x) :> = — 



1 V 



47T/ 



^ 16A 2 K 2 e -2 K (u+B) 

48 ~~ (v - B + Ae~< u + B )) A 



Similarly, the expectation value of the flux is 



1 



2 ^4 



471-/ 48' 



u — ► oo. 



(36) 



(:T ri (x) :) = -(:T H (x) :) = — 



48vr 



The squared flux is related to the mean flux by 

(: T? t (x) :) = 3(: T rt :) 2 , 

and the relative deviation is 



(: AT rt :) _ T? t (x) :) - (: T rt :)2 



V2. 



(37) 



(38) 



(39) 



(: T rt :) (: T rt :) 

The fractional flux fluctuations are thus of order unity. 

Correlation function 

The function {T r t(x)T r t(x')) is finite except in the short distance limit x' — x — > 0. On the 
other hand, the normal-ordered function (: T r t(x)T r t(x') :) has a finite value in this limit. 
Here we restrict our discussion to the latter function. We define a normalized correlation 
function as (note that this is distinct from the function C defined in Eq. (|7j) ) 



£(At) 



(:T H (t)T rt (t f ) :) - (: T rt (t) :) (: T rt {t') :) 



(40) 



(:T rt (t):)(:T rt (t>):) 

where At = \t — t'\ and the spatial points are coincident. For the specific trajectory of 
Eq. pi), £(At) is 



f (At) = 288 



(1 



-K,At 



-nAt\2 



1 



( K At) 2 



(41) 
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Figure 2: The correlation function for thermal radiation in two dimensional spacetime. 
which is plotted in Fig. 2. Note that £(At) is finite for all At > 0. 

e(Te) = ^(0), (42) 

and is approximately 

r c «-. (43) 

K 

3.3 Two dimensional black hole 

Fulling and Davies have shown that the mirror trajectory of Eq. ( |34|) produces the same 
quantum state in the asymptotic region as does a two dimensional evaporating black hole 
of mass M if ac = 1/4M. Thus the flux and squared flux for the 2-D black hole are 



(: T rt (x) :) = 3_, (44) 



and 



(: Tl (x) :) = 3(: T rt {x) :) 2 . (45) 
The correlation time r c becomes 

r c = 11M. (46) 

Thus the Hawking flux undergoes large fluctuations, varying by a factor of order unity on 
a time scale of order 11 M. 
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3.4 Four dimensional black hole 



In four dimensions, the treatment of black hole evaporation becomes more complicated than 
in two dimensions. This is due to the angular degrees of freedom and the resultant potential 
barrier around the black hole. Ingoing and outgoing waves experience scattering off of this 
potential barrier. We will consider the case of a nonrotating, uncharged (Schwarzschild) 



black hole, and will follow the treatment of DeWitt [12]. The mode functions are of the 
form 

u = - l -= R l (u\r)Y lm (e, 0) e~ iw \ (47) 



27rrV2o; 

where the Yi m (6, 4>) are the usual spherical harmonics. Vector signs will be used to indicate 
the two independent modes which have the asymptotic forms 

/ , s e + Anujje ,r — > — oo , ln , 

AHe J * ,r^oo (48) 



and 



<- , , . Bi{uj)e- luir ,r* -> -oo ri . 

e -^V% { ^- r-^oo • (49) 

where r* is the usual tortoise coordinate 

r* = r + 2Mln^-l). (50) 

The transmission coefficients B[ and reflection coefficients A\ satisfy the relations 

\%{uS)\ = \%{^% (51) 

= \-\%{u)\ 2 = \B l {uS)f, (52) 

and 

X{u)B x (u) = -Bf%(oj). (53) 

The components of the energy-momentum tensor in the Unruh vacuum state near future 
null infinity are of the form (See Ref. 12 for details.) 

(<t>,»<j>,v) ~ Y] / + coth(47rMw)l?^"r?* I/ ]dw. (54) 

l,m J ° 

The asymptotic form of the mode functions are 

1? = - ^ l^e^'^ (55) 



and 

= -J-=Y lme -**( e -iur' + ^( W ) e ^). (56) 
27rrv2cj 
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The derivatives of these mode functions become 

■j\Ylm\ 

8n 2 r 2 



frAS = ^^(l-\%\ 2 + %*e~ 2i " r -%e^ r ), (57) 



H^t = ^^(l-fcf-Xe-^ + Xe 2 ^), (58) 



8vr 2 r 2 

lT jt tr* t = ^4(l + \%\ 2 + %*e- 2 ^' + %e 2 ^), (59) 



^ 12 

87T 2 r 2 



w, r V; = ' l m l (1 + \Ai\ 2 + ^ e~ 2lujr + Aie 2wr ), (60) 



— > -** -u\Bi\ 2 \Y lm \ 2 

U ^ U ,t = ^2-2 ' ( 61 ) 

= ^2~2 ' ( 62 ) 

".^.r = 87r 2 r 2 » (63) 

_> _>* w|5ji| 2 |yj m | 2 

= 8n 2 r 2 ' ( 64 ) 



and 



Substitution of these relation into Eq. fl54| ) and use of the summation formula 

£lW = ^, (65) 

m 

yields 



Lm 



/•oo 

22 / [^",r^*"*t + coth(47r M oj)Tt ir Tt* t ]duj 
-sing ^ |vr 2 f°° \B t \ 2 u , 



Lm 



+f^S E( 2 * + 1) /f V*" - V*")** (66) 

The first term on the right hand side is (T r t). Similar calculations give us these derivatives 
in terms of the mean flux (T rt ) as 



\r<f>,r) 



-(T rt )+I 2 , (67) 

-(T tr )-h, (68) 
(Trt) + Jb _j (69) 



and 

where the integrals Iq, I\ and I2 are 

sin# 



(T rt )+I + h. (70) 
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and 

sin 

h = -lT7T 



sin r°° 

— ^V(2Z + 1) / uA(u) cos(o» + 2ru)dw, (72) 

WIT 6 i Jo 

roc 

J2( 21 + 1) J uA{u) siii(<5(w) + 2ruj)dw. (73) 



We can identity Jo as the Minkowski divergence, and symmetrization removes the pure 
imaginary term I<i. Discarding these two terms yields 

i«: 4>,r4>,t ■■) + (■■ 4>,t<t>,r :)) = (:T rt :) = (:T tr :), (74) 
(: Mr ■) = -(:T rt :)-I 1: (75) 

and 

(:^, t :)=-(:T rt :)+/!. (76) 

Using the relations 

A{u) -> A , uj^ 0, 
A(uj) — > 0, w — > oo, 

and 

5(w) — > <5o, u; — ► oo, 
the integral ii at large distance, r —>■ oo, becomes 



1 Z" 00 /x . 

xA — cos 



<5 ( - ] + 2x 

r 



rOO Y f 

/ ujAUjj) cos(<5(a>) + 2ruj)duj = 

Jo r 2 Jo 

~ — =- / cos(2x + oo )xax oc 



du> 



(77) 



Thus the term I\ is much smaller than the mean flux 

T 



which is a nonzero constant at large distance, and hence I\ can be ignored. The squared 
flux becomes 

( : T rt ■) = (■ T rt :) 2 + (: <\>At '■)('■ $,t<$>,r ■) + (: 4>,r4>,r ■}(■ 4>,t<t>,t ■} 
= 2<: T rt ■? + ((■■ T rt :) - T rt :) + h) 

~ 3<T rt > 2 . (79) 

We thus get the same relation between the squared flux and the mean flux as in the case 
of a two dimensional evaporating black hole. 
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Next we will discuss the normal-ordered correlation function and ignore the contribution 
from I\. From Eq. (|54|), we have 



and 



(■ <t>,r(x)4>,t(x') : 
(: 4> t t(x)4>, r (x') 



Re 



Re 



Re 



16vr 



16vr 



87rA/o; 



00 IBiI'W 



1 r 

^E( 2Z + 1 ) / 

T 2 r 2 Jo 

1 POO I D 

1 /" 



00 iSjPwe^* 



1 



-do; 



duj 



B t \ 2 u;e iujAt 



16vr 2 



Re 



16vr 2 



e 8ivMu> _ J 



?8ttMu) _ I 



cLu 



(80) 
(81) 
(82) 

(83) 



where At = t' — t. The correlation function defined in Eq. (f40|) now becomes 



£(At) 



2(ite [Ez(2/ + 1)/ C 



-1 



-do; 



(84) 



(ite [E,(2l + 1) /o 00 j£r^ 
For our purposes, the transmission coefficient may be approximated by a step function 

Bi{p) « e(V27Ma; - Z). (85) 



In this approximation, modes with energies below the peak of the angular momentum 
barrier are assumed to be perfectly reflected, and those with energies above the peak are 
completely transmitted. This is a reasonably good approximation, as may be seen by 
examining Figure 1 in Ref. [13], where numerical results for the transmission and reflection 
coefficients are given. The summation on I yields 



^(2Z + l)\Bi\ 2 = ^2(21 + l)9(\/27 Muj - 1) = 27MV + 6V3Mu + 1 
1 1 



(86) 



and the numerator of Eq. ( [B4] ) becomes 

(poo I R,|2 iuAt 



dw 



Re 



(27M 2 uj 2 + 6y/3Muj + l)u;e ipAt 



27 



3vr) 4 M 2 
1 



+- 



Re 



e 8irMuj _ ^ 

da; I + 



00 x 3 e ibx 



6^/3 



e 



' - 1 

ibx 



3vr) 3 M 2 



da; 



00 x^gibx 



dx 



e J 



-dx 



(8ir) 2 M 2 ~~~ \J e x -l 
where x = 8ttMuj and b = At/8nM. The correlation function is 

2S(At) 2 



(87) 



£(At) 



S(0)= 



12 



2 




At/S-rrM 



Figure 3: The correlation function for the radiation from a black hole in four dimensional 
spacetime. 



and is plotted in Fig. 3. 

The correlation time r c is around 0.3(87rM) ~ 8M . As in the two dimensional case, 
the four dimensional Hawking radiation undergoes large flux fluctuations on a time scale 
of about 8M. 



3.5 Flux Fluctuations as Thermal Fluctuations 

It is reasonable to expect that the flux fluctuations computed in the previous subsection 
can be interpreted as ordinary thermal fluctuations. Thermal fluctuations of the energy in 
the canonical ensemble are described by the relation 



AE _ T^C^ 
E E 

where E is the mean energy at temperature T for a system with heat capacity Cy. In the 
case of thermal radiation, E oc T 4 , so Cy = AE/T, and hence 



AE k B T , . 

Note that E/iksT) is a measure of the mean number of photons in the thermal radiation, 
so the above result is the familiar 1/y/N statistical fluctuation. 

Let us take E to be the energy emitted by a black hole in one correlation time, r c ~ 8M, 
and the power emitted to be that calculated by Page 14] for photon emission from a 
Schwarzscild black hole: 

P = 3.4 x Kr 5 M" 2 . (91) 
This leads to a rather small number of photons emitted per correlation time, 

TP 

3 x 1(T 3 (92) 
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and rather large fractional energy fluctuations 



AE 

— « 23. (93) 

This estimate is somewhat larger than the result obtained from the normal ordered 
squared flux in the previous subsection. However, it is a very rough estimate which depends 
upon our choice for the energy E and the collecting time. If we had chosen to integrate 
the flux for a time longer than r c , the statistical fluctuations would be somewhat reduced. 
Also, the spectrum of particles emitted by a black hole is not exactly Planckian, but has 
been filtered by the angular momentum barrier. 



4 The Physics of the Cross Term 

We now turn to examining the cross term between the vacuum fluctuations and the finite, 
state dependent parts. Recall that this contribution is singular in the coincidence limit, 
and hence does not lead to a well defined local definition of the squared flux. If it is not to 
be subtracted by some renormalization method, then it can only be given physical meaning 
by dealing with time or space averages. 

4.1 Switching Functions 

One possibility is to suppose that we operationally measure the flux with a model detector 
which has a finite response time. Suppose that the response of our detector is described by 
a Lorentzian function with characteristic width r, 

m = ; (t _ to ; 2+T2 . 

The averaged squared flux becomes 

/oo roo 
/ f{t)f{t'){T rt {x)T rt {x'))dtdt'. (95) 
-OO J —CO 

We will examine the case of the thermal flux from a mirror or black hole in two dimensions, 
and assume that the sampling time is short compared to the correlation time, r <C r c . 
In this case, the correlation functions are approximately constant, and the average of the 
normal ordered term is 

/oo /*oo 
/ f(t)f(t'){: T rt {x)T rt (x') :}dtdt' « (: T rt {x)T rt {x') :). (96) 
-OO J —CO 

The average of the cross term can be written as 

/oo rco 
/ f(t)f(t') [(: </>(x), r <f>(x'), r , ■)(<t>{x) :t <t>{x'), t ,) M + 
-CO J —CO 

/CO 1 

' (T rt ), (97) 



4ttt 2 
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where we have used Eq. (B8) in Appendix B. 



For comparison, it is of interest to give the average of the vacuum term 

/*oo /*oo 

A) = / / /(t)/(i / )(^(^),r0(x , ),r')A/(0(^),t0(x , ), t ')M^ / 



OO J — oo 

■ poo i i 



mf ^4^W d ' dt ' = M^' (98) 
using Eq. ( gig) . 

In the case of a two-dimensional black hole, where (T rt ) = l/(7687rM 2 ) and (: T 2 t :) = 
3(T rt ) 2 , we see that 3> A\ 3> ^ when T 3> M. If we were to let T « M, our assumption 
that the correlation functions are constant would no longer be exact. Nonetheless, this 
calculation should give a resonable order of magnitude estimate, and predicts that all three 
terms are of the same order of magnitude. 

4.2 A Mirror as a Flux Detector 

Here we will examine a model in which the flux is measured by the force which it exerts on 
a reflecting or absorbing surface. All of our discussion will be in one spatial dimension, so 
the force on a partially reflecting surface is 

F = r : T rt : , (99) 

where < r < 2 and \r is the fraction of the radiation which is reflected. Consider a mirror 
with mass m which starts at rest at time t = . The mean velocity and mean squared 
velocity at t = r are 

(v) = — f T (F)dt (100) 



m Jo 

and 



(v 2 ) = \\ I (F(t)F(t'))dtdt'. (101) 



m 2 Jo Jo 
The force two-point function is 

(F(t)F(t')} = r 2 (:T rt (t)::T rt (t'):} 

= r 2 (: T rt (t)T rt (t') :) + r 2 (T rt T rt ) cross + r 2 (T rt T rt ) M . 

(102) 

We are interested in the difference between the fluctuations in a given state and those in 
the Minkowski vacuum, and so drop the vacuum term. We then define the fluctuations by 
subtracting out the square of the mean value: 

(AT 2 ) = (T 2 ) - (T) 2 . (103) 

The normal-ordered flux fluctuation is given by 

(AT 2 ) wo = (:T rt (t)T rt (t>):)-(T rt (t))(T rt (t>)) 

= (: 0, r (t)£ r (f) :)(: <M^(0 :) + <: 0,r(*)^t(f) :)(: MWA*) : >> 

(104) 
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and the cross-term contribution is 

(AT 2 ) cross = (: <M*)<MO WAWMm + <: MWttf) 0<<M*)<W))m. (105) 
The velocity fluctuation can be written as 

(At, 2 ) = (v 2 )-(v) 2 

= (Av 2 ) NO + (v 2 ) 



cross 

2 r Jl 



m? 



f (AT 2 t ) NO dtdt' + ^ / (AT 2 ) cross dtdt'. (106) 

J 7Tb J 



4.2.1 Coherent state 



A coherent state describes a classical field excitation and is hence a useful model to reveal 
the effects of the cross term. Consider a single-mode coherent state \z) for a mode with 
frequency ujq 

a w \z) = 5 UUo z\z). (107) 
The free quantum field expanded in normal modes is 

where the mode function for a standing wave in a box of length L is 

' e-^sin^x). (109) 



Assume that the mirror remains approximately stationary, as will be the case when its mass 
is large, and set x = x' . Further, let z = Re 1 and find 

(z\ : 4>^{x)4>y{x') : \z) = ^\z\ 2 cos 2 {u: r) _ ^ + cog(2A _ ^ + ^ (nQ) 

and 

(z\ : 4>, t {x)(j>, t ,(x') : \z) = ^ Sin ^° rj [ C o S (LO G {t-t'))-co S (2A-2uj a t+uj a {t+t'))). (Ill) 

Li 

We also have that 

(((>(x),t<f>(x'),t>)M = {(f>( x ),r<t>( x '),r')M = ~ 2 Tr(t-t') 2 ' 

For the coherent state, the only fluctuations come from the cross term because 

(z\ : T^ v T pa : \z) = (z\ : : \z)(z\ : T pa : \z). (113) 
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The velocity fluctuation is then 

(At. 2 ) = ^ I {AT 2 t ) cross dtdt' 



-uj \z\ 2 r 2 /' 2 



irLm 2 



{cos (LOor)[cos(ioo(t — t')) + cos(2A — uo(t + t'))] 



+ sin 2 (cj r)[cos(u;o(t - t')) - cos(2A - 2u t + u (t + t' ))]}-, — ^—^dtdt'. 

{t - t'y 

(114) 

This integral is poorly defined due to the singularity of the integrand at t = t' . A 
possible resolution of this difficulty is to employ a trick which has been used by various 
authors under the rubrics of "generalized principle value" ][l5|] or "differential regularization" 
1 16]. In any case, it involves writing the singular factor as a derivative of a less singular 
function, and then integrating by parts. We will also assume that the flux is adiabatically 
switched on in the past and off again in the future, so that any surface terms vanish. Thus 
we may use relations such as 

/ fMr^iydtdt' = l -j [d t d t ,f(t,t')} ln((t-t') 2 )dtdt>. (115) 

Let r = , x = toot, and x' = LO$t' . The velocity fluctuation becomes 

(Av 2 ) = - r ( ( [cos(u t-u t')+cos(2A-uJot-u;ot')}ln[LOo(t-t')] 2 dtdt' (116) 

l-nLra 1 Jo Jo 

The integral may now be written in terms of the variables U = x — x' and V = x + x' , using 
the identity 

1 / f° fU+2r rr f2r-U \ 

dtdt' = - dU dV+ dU dV , (117) 



o Jo 2 \J-t J-u Jo Ju 

and evaluated in terms of sine and cosine integral functions. We are primarily interested 
in the asymptotic form for large r, which is 

6Lm z 

This result shows that the cross term leads to a contribution to the mean squared veloc- 
ity of the mirror which grows linearly in time. This is the characteristic time dependence 
of a random walk process. It is useful to compare the fluctuations with the mean velocity 

(v) = -f T (T rt )dt 
m Jo 

I 1 2 

= ^-^sin(2a;oa;)[cos(2A) - cos(2A - 2u t)]. (119) 

The mean velocity happens to vanish at the special point x = at which we evaluated 
(Av 2 ). However, at a more general point it is of order 

I 1 2 

(: „ :) - LL_, (120) 
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and the fractional velocity fluctuations becomes of order 



(121) 



Although this quantity grows in time, it is also inversely proportional to the amplitude \z\. 
Thus for a nearly classical state (\z\ 3> 1), it can remain small for a very long time. 

4.2.2 Thermal state created by moving mirror 

Consider the 2-D moving mirror with an arbitrary trajectory. These flux fluctuation, with- 
out the vacuum term, can be written as 



(AT r 2 



47T 



rt I total 



1 



<Ar 2 > no + (AT 2 



where 



(AT 2 



rt/no 



(p(u) — p(u')) 2 
p'{u') 



•(- 



rt I cross 

gV) 

[v — p(u')) 2 
p'(u) 



) 2 -( 



p'(u) 



{v-p(u)) 2 > (t-f) 4 
p'(u)p'(u' 



(v — p{u')) 2 
2p'{u)p'{u') 
(p(u)-p{u')) 2 {t-t') 2 



' (v — p(u)) r " 



+ 



(p(u) — p(u')) 2 



r + 



(t - tr 

(123) 



and 



(AT 2 ) 



p'{u)p'(u') 1 

(p(u) -p(u>)) 2 ~ (t - t'f 



(t - 1>) 



A2- 



(124) 



The total velocity fluctuation is 



total 



— / / {AT? t ) total dtdt' 



Jo Jo 
2r 2 

(4ir) 2 m 2 Jo Jo 



p'(u)p'(u') 



i 2 



(p(u) — p(u')) 2 



p'(u) 


2 


r P '(u') - 


2 - 1 1 


(v' — p(u)) 2 _ 




. ( v — p{u')) 2 _ 


(t-trj 



dtdt'. 



(125) 



As before, we suppose that (AT 2 t ) to tai is multiplied by a swithching function which vanishes 
in the past and in the future. Then the (t — t')~ 4 term gives no contribution after an inte- 
gration by parts. Note that this term is of the same form as the pure vacuum contribution, 
so our results will not depend upon whether the vacuum part was subtracted beforehand 
or not. 

Consider the trajectory which produces a thermal spectrum, p(u) = B — A e ~ K ( u + B ) anc l 
again assume that the detector remains at a fixed location, which we take to be x = 0. 
Then u = v = t and v! = v' = if . The integral of two middle terms in Eq. ( |125| ) may be 
shown to approach a constant in the limit of large sampling time r: 



o Jo 



p'{u) 



i 2 



(v' — p(u)) 2 



+ 



p>(u>) 



{v — p{u')Y A 



dtdt' 



constant. 



(126) 
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As we will see, this is small compared to the leading term, which grows linearly in r. The 
total velocity fluctuation is now 



{Av 2 ) tota i ~ 



2r 2 



(47r) 2 m 2 Jo Jo 



p'(u)p'(u') 



n A r 2 



87r 2 m 2 7o (l-e K (*-*')) 



(p(u) -p{u')) 2 

2n(t-t') 

dtdt' = 



dtdt 



k a r 2 
4iT 2 m 2 



2kU 



(t-U) 



(l_ e «tf)4 



dU, (127) 



where in the last step we have used the change of variables given in Eq. ( |117| ). At late 
times, we have 



(A*, 2 ) 



4 2 
TK r 



total 



„2kU 



47r 2 m 2 Jo {l-e KjJ ) A 



dU 



S 2 

t k r 
2TT 2 m 2 



n 2x 



(1 



x\4 



dx , 



(128) 



where x = kU. If we ignore the singularity in the integrand, this integral may be evaluated 
directly: 



•2x 



1 



x\4 



dx 



3e a 



1 



6(1 



oX\3 



1 

6 



A more rigorous approach is to use the relation 

1 1 d A 



Inx 



x 4 12 dx A 

integrate by parts, and then evaluate the resulting integral numerically: 



1 



2.r 



-oo x 4 (1 - e 
In either case, the result is 



dx = 



12 J- 



(Av 2 (T)) total 



dx 4 



x A e 2x 



(1 



oX\4 



dx w - . 



k 3 r 2 



(129) 



(130) 



(131) 



(132) 



487r 2 m 2 

As in the case of the coherent state discussed in the previous subsection, the mean squared 
velocity fluctuations grow linearly in time. 

We now wish to determine the relative contributions of the normal-ordered and cross 
terms. A calculation analogous to that performed for (Av 2 ) total reveals that the normal- 
ordered contribution is also linearly growing in time: 



(A« 2 ) r 



tk 3 r 2 
2TT 2 m 2 



1 

~ A + 



-,2x 



(1 



oX\4 



X 2 1 



oX\2 



dx . 



(133) 



In this case, the integrand is finite from the beginning, so no integration by parts is needed. 
The integral may be evaluated numerically to yield 



(A^ 2 ) nc 



k 3 r 2 T 
192vr 2 m 2 



(1.00). 



(134) 
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The cross term contribution may be obtained as the difference of Eqs. ( |132j ) and (|134| ), 
but it is useful as a check to compute it independently. If we follow the procedure used to 
find the asymptotic form of (Av 2 )totah including an integration by parts, the result is 

Again, the integral may be evaluated numerically with the result 

3 2 

{Av 2 ) cross ~ " I T 2 (1.00) « 3{AT 2 t ) no . (136) 

Thus to the accuracy of the numerical calculations, the three independently computed 
pieces do indeed satisfy 

(Av 2 ) total = (Av 2 ) no + (Av 2 ) cross . (137) 

The normal-ordered term contributes 25% of the total effect, as compared to 75% from the 
cross term. 



6 e x {e 2x + 4e x + 1) 



[1 



dx . 



(135) 



4.2.3 Mass Fluctuations of Two Dimensional Black Holes 

We may use the above results to discuss the fluctuations in the mass of evaporating black 
holes in two dimensions. Define the mass operator by 

M(T) = M - / T T rt {t)dt , (138) 
Jo 

where Mq is the initial mass at time t = 0. The mean mass decreases according to the 
semiclassical theory of gravity as 

(M) = M - [ T (T rt (t))dt, (139) 
Jo 

However, the squared mass will undergo fluctuations: 

(M 2 ) - (M) 2 = [ T f T [(T rt (t)T rt (t')) - (T rt (t))(T rt (t'))]dtdt' . (140) 
Jo Jo 



Equations ( 125 ) and (|132j) may be used to show that 

(M 2 ) - (M) 2 ~ ^ « T -^, (141) 

V 1 V 1 48vr 2 1152M 3 ' V ; 

where in the last step we assumed that n = l/(4Mo), which is a good approximation in the 
early stages of evaporation. We may estimate the evaporation time T evap of a black hole by 
setting 

Mo , , 

Tevap « /rr = 768vrM 3 . (142) 



(Th(0)) 
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If we set r = r evap in Eq. (141), the result is 



(M 2 ) - (M) 2 ps — . (143) 

3"7T 



Recall that we are working in Planck units, so the right hand side of Eq. ( |143|) represents 
a mass fluctuation of the order of the Planck mass. Even though this effect is quite small 
for macroscopic black holes, the cross term plays a significant role here. If one were to 



normal order the product of flux operators above, then the right hand sides of Eqs. (141) 
and ( |143[) will be decreased by factor of 1/4. This leads to a thought experiment in which 
one could use evaporating black holes to test the reality of the cross term. One prepares 
several black holes with the same initial mass, and then measures the masses at some later 
time. If the mass fluctuation grows linearly in time in accordance with Eq. ( |141| ) (or its 
four dimensional analog), then one would have measued the effect of the cross term. 



5 Conclusion and discussion 

We have seen that the flux of radiation from an evaporating black hole undergoes large 
fluctuations on short times scales. One approach to the subject of flux flucuations involves 
the use of normal ordered products of stress tensor operators. In this approach the squared 
flux is a finite, local quantity and the Hawking flux undergoes fluctuations of order one 
on times scales of order M, the black hole's mass. These fluctuations can be viewed as 
esentially statistical fluctuations due to the small mean number of particle emitted by the 
black hole on this time scale. 

However, the subject of stress tensor fluctuations can be a subtle one, and there is 
another approach in which one retains the state dependent cross term in the product of 
stress tensor operators. This term is divergent in the limit that both operators are evaluated 
at the same point. Consequently, if it is present one cannot give a meaning to the local 
squared flux. It is, however, still possible to define time integrals of a product of fluxes. 
These time integrals may be used to show that, at least in a two dimensional model, the 
mean square mass of a black hole differs from the square of the mean mass by an amount 
which grows linearly in time. Furthermore, the rate of growth is four times larger when 
the cross term is retained as compared to the normal ordering approach. Thus in principle, 
the two approaches have different observational consequences. Similarly, they give different 
predictions for the velocity fluctuations of a material body, such as the mirror discussed in 



Sect. 4.2 



In either approach, we are dealing with thermal radiation (or filtered radiation) in the 
asymptotic region far from the black hole. Thus the ambiguity in how to treat the fluc- 
tuations is not confined to the specific case of a black hole, but is present in a general 
quantum state. Because we are working in the asymptotic region, we cannot directly ad- 



dress the issue of horizon fluctuations caused by quantum stress tensor fluctuations [17]. 
Horizon fluctuations must be far below the Planck scale in order that Hawking's semiclas- 
sical derivation Q of black hole evaporation hold. Estimates of the scale of the scale of the 
horizon fluctuations due to quantization of the gravitational field ("active" as opposed to 
"passive" fluctuations) indicate that Hawking's derivation does indeed hold for black holes 
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above the Planck mass |Tq] . It is thus of interest to calculate more carefully the scale of 
"passive" fluctuations due to stress tensor fluctuations. 

Acknowledgement: This work was supported in part by the National Science Foun- 
dation (Grant No. PHY-9800965). 

Appendix A 

Assume there is a quantum state \ip) a which can decompose the field operator into positive 
and negative frequency parts = + + 0~, with (j) + \ip) a = 0. By using Wick's theorem, 
the four-point function can be expressed as 

01020304 = N a (<f>l 020304) + Aa(0102)(0304)a + N a (<j>l 03) (0204)a + N a (01 4 ) (0203)a 
+ A r a(020 3 )(0104>a + iVa(0204) (010 3 )a + N a (0 3 4 ) (01 2 ) a + (010 2 )a(0304>a 
+ (0103)a(0204)a + (0104)a(0203)a , (Al) 

where N a means normal ordering with respect to the state \vp) a and () a means the expecta- 
tion value in this state. If we take the expectation value of the above relation in the state 
\ip) a , the result is 

(01020304 )o = (0102)a(0304)a + (0103)a(0204)a + (0104)a(0203)a- (A2) 

For the particular case of the Minkowski vacuum, 

01020304 = : 01020304 ■ + ■ 0102 = (0304)*/+ ■ 0103 ■ (0204>Af 

+ : 0104 : (0203)a/+ : 0203 : (0i04)a/+ : 0204 : (0103)m 

+ : 0304 : (0102)*/ + (0102)m(0304)m 

+ (0103)Af(0204)Af + (0104)A/(0203)A/ , (A3) 

where :: means normal order to the Minkowski vacuum state, and {)m means the expecta- 
tion value in the Minkowski vacuum. The expectation value of the above equation in this 
state is 

(01020304)Af = (0102)A/(0304)Af + (01 03) M (0204) M + (0104) Af (0203>Af • (A4) 

By using the expressions 

0102 =: 0102 : +(0i02)a/ (A5) 

and 

(0102)a = (: 0102 :)a + (0102>Af , (A6) 

we get 

(: 0102 :)o = (0102)a - (0102)Af (A7) 

and 

(: 01020304 :)a = (: 0102 :)a(: 0304 :)a + ( : 0103 : )a(: 0204 -)a + (- 0104 -)a(- 0203 -)a ■ (A8) 
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Appendix B 

Bl Evaluation of J f(t)f(t' )-^ F ^dtdt / 

The sampling function and its derivative are respectively 

r 1 

/(*) = z; 



vr(t-t ) 2 +T 2 
and 

fit) = ^2 

/U vr ((t-t ) 2 +T 2 ) 2 ' 
The integral after intergation by parts yields 



where 



f{t)f{t')j^dtdt' = J J f{t)f'{t')Ht-t')dtdt' 

/roo 
f'(t')dt> f'(t)ln(t-t')dt 
J —CO 

= I f'(t')A(t')dt', 



((t-t ) 2 + T 2 ) 2 

This integral contains a second order pole and can be done by residues. Let x = t 
the integral becomes 

A(t') = — j ^ (g2 + ^ 2)2 ln(x + 1 - t')dx 

= f ~nr t — tv? ln ( z + *o - * )dz 

= =T (" w) f ?TF 1,1(3 + t °- 

-(27Tt) 



7T 





■ln(z + t -O" 


(dr 2 ) 


Z + IT 



b t 
+ i- 



T 2 +b 2 V+6 2 ' 

where R = \/t 2 + b 2 and b = t' — to • The double integral becomes 

-2r 7r[4r 2 -(t -t ) 2 ] 



vr 2r[4r 2 + (t - t ) 2 ] 2 
We keep only the real part and write 

-2r vr[4r 2 - (t - t ) 2 ] 



/'(tO/m^^Oldt' 



(t-t') 2 vr 2r[4r2 + (t -t ) 2 ] 2 ' 

In the coincidence limit to t' Q , the integral is 
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B2 Evaluation of / f(t)f(t')j^dtdt' 
The second derivative of sampling function is 

1 



/"(*) = -- 

7T 



Integrating by parts yields 



4(t - t f 



((t-t ) 2 + r 2 ) 2 ((t-t )2 + r 2^3 



f{t)f{t')-^y 4 dtdt' = -1 I I f{t)"f{t')"Ht - t')dtdt'. 



(B9) 



(BIO) 



A similar calculation as in Bl yields 

POO 

B{t') = I f"{t)\n(t-t')dt 

\xiz 



2r 

7T 

2a 

TT 



_d_ 

dr 2 I (z-f- t ) 2 + r 2 
B 2 



dz + 2 



( d \ 2 r (z + f-t f 



\dT 2 J J (Z + f - to) 2 + T 2 



In zdz 



d „ / 3 \ 2 
— ? 5i + 2 — - 



<9r 2 



\dT 2 J 



where 



and 



Bi = — In zq 

T 



B 2 = -7rrlnz , 

where zq = \Jt 2 + b 2 and b = t' — to- The integral B(t') becomes 



B(t>) = J f"{t) Ht - t')dt = --^ + j- 



2r 2 



(r 2 + 6 2 ) 2 



(Bll) 

(B12) 
(B13) 

(B14) 



Substituting this result into the original double integral yields 



f(t)f(t')- 



(t - tr 

-t r°° l 

3^ J-JW -tn) : 



dtdt' 



4(f - t'o) 2 . . 1 



r)(: 



2r 2 



, 3 ((tf - t' ) 2 + r 2 ) 2 ((t'-^) 2 + r 2 ) 3M r 2 + 6 2 (r 2 + 6 2 ) 2 
16r 2 -24r 2 (t -t , ) 2 + (to-^) 4 



)dt' 



(Ar 2 + (to - t' ) 2 r 
In the limit to — t' — > 0, the integral becomes 



(B15) 



(B16) 
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